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I. INTRODUCTION AND OVERVIEW 

System identification technology has been used 
successfully for many vehicles. Because of their 
large number of degrees of freedom and complex aero- 
dynamic interactions, the rotorcraft have always 
presented a special challenge to system identification 
methods. A completely integrated methodology has been 
developed under this NASA contract to solve this 
difficult problem. This methodology has also been 
translated into a user oriented series of computer 
programs. This volume provides basic guidelines for 
efficient and effective use of one of these computer 
programs . 

Figure 1 shows a schematic flowchart of the overall 
data processing technique for rotorcraft. The first 
step in this procedure is state estimation and instrument 
calibration. This is implemented by the computer program 
DEKFIS (for Discrete Extended Kalman Filter and Smoother 
which implements an extended Kalman filter/smoother using 
the Friedland-Duffy formulation. Instrument biases and 
scale factors are estimated at this stage together with 
any state which is not measured directly. The second 
step involves. estimation of the mathematical model of 
various forces, moments and interchanges. This is 
implemented in OSR (Optimal Subset Regression) computer 
program which uses a regression technique. Accurate 
estimates of parameters are obtained in the final step. 

One of two computer programs is used for this purpose. 
SCIDNT implements the maximum likelihood method for linear 
systems and NLSCIDNT extends the method to nonlinear 
rotorcraft models. 


The contract research effort which led to the results in this report was 
financially supported by the Structures Laboratory, USARTL (AVRADCOM), 
NASA Langley Research Center and NASA Ames Research Center. 


1 




Figure 1 Integrated Rqtorcraft System 
Identification Procedure 


Accuracy of parameter estimates may be improved 
by using flight test inputs based on the input design 
program, INDES. 

This userts manual describes the OSR 
computer program. The details of the theory and the 
particular implementation used are given in the final 
report . * 


* Hall, W.E., Gupta, N.K., Hansen, R. and Bohn, J., 
’’State Estimation and Parameter Identification 
for Rotorcraft,” final report on contract NASI-14549, 
May 1978. 
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II. BACKGROUND 


The computer program, OSR (Optimal Subset Regression) 
estimates models for rotorcraft body and rotor force and 
moment coefficients. The technique used is based on the 
subset regression algorithm, described extensively in 
Draper and Smith [1]. Given time histories of aerodynamic 
coefficients (e.g. C^, C^, etc.), aerodynamic variables 

(e.g. forward speed, angle-of-attack, etc.) and control inputs, 
the program computes correlation between various time histories 
and the model structure determination is based on these cor- 
relations. The procedure may be described as follows. Consider 
a single force or moment coefficient y and assume that it 
may be a non-linear function of several aerodyanmic variables 

or control inputs, x' (called independent variables). A 

number of functions of x" are formed to describe all 
possible relations between y and x". Let x be a 
collection of such functions; y and x are assumed to 
be related by the equation: 

y - x^e . e (1) 

where 0 is anMxl vector of unknown parameters and e is 
noise. For sake of simplicity e is assumed white and 
gaussian. If there are N measurements. Equation (1) 
may be written as: 

Y = X0 + e 

where y is an Nxl vector and x is an NxM matrix. The 
program then performs a series of hypotehsis tests based 
on the F-ratio criterion to determine which 9's are zero. 

The X terms corresponding to non-zero 0’s are the model 
forms for the corresponding aerodynamic force or moment 
coefficient. 

The program has the following features : 



(1) The program is set up for rotorcraft models. 

(2) The coefficients, x, are specified in a separate 
subroutine. Therefore, the form and number of 
functional forms of basic variables may be 
changed easily. 

(3) Each force or moment coefficient is treated 
separately, allowing maximum flexibility. 

(4) lipto 80 terms may be specified in each model. 

(5) Any term may be forced in the equation. Also 
terms may be forced out of the equation. 

C6) Time histories of measured and estimated variables 
(based on the model structure) may be printed 
out if desired. 
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III. PROGRAM STRUCTURE 


This section gives a brief discussion of the major 
subroutines and functional blocks of OSR. Figure 3.1 
shows the basic program structure. 

The main routine, MAIN reads the data sets up the 
problem and does most of the computation. Most of the other 
subroutines are special purpose subprograms. From user 
viewpoint only the TRANSG subroutine is of importance. 

The TRANSG routine uses an initial set of variables 

to create new sets of variables. For example given a 

2 

we may create an a time history in this subroutine. Other 
variables may be added as desired. In general this sub- 
routine has to be changed for each run. However, for projects 
in which the model (i.e., TPJ*JISG) remains unchanged, the 
appropriate TRANSG subroutine can be incorporated into OSR 
at the beginning of the project and remain unchanged through 
the end of the project. 

Other parts of Figure 3.1 are self descriptive. 
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MAIN 


Read input 
1 


Create new variables 


RDBL2 

Read and write labels for printout 
TRANSG 

Compute new variables based on input data 


Compute 'covariance matrix 

Compute^correlation matrix 

Enter a^variable (sTEPRG 

I I Determine variable to 

^be included 


Compute regression equation coefficients, 
standard error and F-ratio-to-remove for 
for variables in the equation 

Compute ^partial correlation 
coefficient, tolerance and 
F-ratio for variables out 
of the equation 

Print regression table jwHICHX 

Check if a variable should be 
removed from the equation 

♦ 

Check if a new variable may be 
entered. If not proceed to next 
step I 

T ( RES IDS 

Compute residuals < Reread mass storage device 

Print risiduals, if desired 


STEP 

Compute parameter 
estimate and update 
correlation matrix 
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IV. INPUTS TO THE PROGRAM 


OSR requires three classes of input. The first type, 
which always consists of cards, defines the structure of 
the time history file, specifies program tolerances and 
constraints, and dotermincs types of output ♦ 

The second type of input, referred to as the time- 
history input or input data file, consists of tabular 
values of control settings, state variable, and other 
relevant model parameter time histories. This data can 
consist of cards, which would then be read with data of 
the first type, or a file on a mass storage device (disk 
or tape) . No choice as to the logical unit number is 
permissible and it must be TAPE4 . The format of the input 
data file is defined in the input of the first type. 

The third type of input ■ is not program data but rather 
a user supplied subroutine. The purpose of this routine 
is to enable the user to generate additional variables 
which are functions of the original variables read from 
the input data file. The structure of the subroutine, 
called TRANSG, is shown in Figure 4.1. The first NOV ele- 
ments of array X contains the original variables read from 
the input data file arranged according to the sorting 
array (see TABLE 4.1, card type 4). 

A sample input deck is shown in Figure 4.2. 
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Table 4.1 


CARO TYPE COLUMNS 

FORMAT 

NAME 

DESCRIPTION 

1 

1-52 

13A4 

ITITLE 

Alphanumeric title for the problem 

2 

1-5 


N 

Number of data points to be used 


6-10 

15 

NSKIP 

Number of data points to be 
skipped on input data file 


n-15 

15 

NOV 

Number of variables from input 
data to be used in regression 
analysis 


16-20 

15 

NVA 

Number of variables added by 
transgeneration (i.e., by 
TRANSG) to be used in regression 


22-25 

A4 

ZERO I 

YES if the y regression equation 
is to pass through the origin, 

NO if not. (The word YES or NO 
must be left justified) 


27-30 

A4 

ST0R9 

YES if the mean and covariance 
matrices are to be written on 
TAPE9, otherwise leave blank 


31-35 

15 

NITEM 

Number of variables read per 
record on input data file 
(NITEM ^ NOV) 


36-40 

15 

NTRANS 

=0 TRANSG hot called and no 
variables are added by trans- 
generation . 7^0 TRANSG is called 

and variables are added by 
trarisgeneration 

3 

Blank card required 


4 

1-80 

4012 

ISORT(I), 
1=1 ,NOV 

Vector to specify the selection 
and order of variables to be used 


in the analysis. The Ith variable 
in the analysis will be the 
ISORT(I) element of the input data 
file 



Table 4.1 (Continued) 


CARD TYPE COLUMNS FORMAT NAME 


DESCRIPTION 


5 

1-3 

13 

NAIT 

Logical unit number for data 
input file. (May be 5 for 
card input or 7 for disk or 
tape only) 


4-6 

13 

NVFC 

Number of variable format cards. 
(This card is the first and 
may be the only one, ten is 
maximum.) 


9-80 

18A4 

RES(I), 

1=1,18 

Format specification for reading 
input data file; if RES(l) is 

6 

1-80 

20A4 

RES(I). 

1=19,... 

blank, file, is assumed to be 
unformatted data 

7 

.1 

A1 

ECHK 

Blank if table information 
follows. * terminates the 
reading of this card type and 
the remainder of this card is 
ignored 


2-3 

12 

K 

Index of the variables cor- 
responding to the first of 





the following tables 


5-76 

9A8 

DUMY 

8-character tables for variables 
K to K+8. NOTE: only six 

characters are printed in most 
of the output 

8 

1-7 . 

1-7 

LB IAS 

=TRUE if any of the variables 
are to be biased 
=FALSE if none are biased and 
no cards of type 9 are read 

9 

1 

A1 

ECHR 

Blank if bias information 


follows. * terminates the 
reading of this card type and 
the remainder of this card 




Table 4,1 (Continued) 


CARD TYPE 

COLUMNS. 

FORMAT 

NAME 

DESCRIPTION 


2-4 

13 

L 

Index of variable to be biased 


5-15 

Fll.O 

BIAS(L) 

Value of the bias 
NOTE: bias will be subtracted 

from variable before trans- 
gression variables are added 

10 


RES 

X(I),I=1, 

NITEM 

If time history data file 
consists of cards read them here 

Subproblem cards; may include as 
desired. 

many subproblem card sets as 

11 

1-48 

12A4 

ITITLE(I), 

1=14,25 

Alphanumeric title information 
identifying subproblem 

12 

1-5 

15 

KDEP 

Index of the dependent variable 


6-10 

15 

MAXSTP 

Maximum number of iterations oer- 
mittedMAXSTP<2*( NOV+NVA) 


13-15 

A3 

RESID 

=YES if residuals are to be 
printed or plotted 


18-20 

A3 

SUMTAB 

=YES if summary table is to be 
printed 


21-30 

FIO.O 

FIN 

F value for inclusion 


31-40 

FIO.O 

FOU 

F value for deletion 


41-50 

FIO.O 

TO 

Tolerance level 


51-60 

FIO.O 

ERTEST 

Regression steps will step when 
r2 >1 .0-ERTEST 


61-65 

L5 

PRTRES 

TRUE if table of residuals is 
to be printed 


66-70 

L5 

PLTRES 

TRUE if plots of residuals are 
desired 

13 

1-80 

8011 

C(I),I=1, 

NOV+NVA 

Flag to indicate status of the 
variables in this subproblem. 
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Table 4.1 (Continued) 


CARD TYPE COLUMNS FORMAT NAME 



DESCRIPTION 


C(I)=1; Ith variable is not to 
be used in any possible re- 
gression equation 

C(I)=2; Ith variable is an 
independent variable for this 
subproblem and may be used in 
regression equation 

C(I)=3,4,. ..,9; Ith variable 
will be forced into the re- 
gression equation. All 
variables with C(I)=9 will be 
forced in first followed by 
variables with C(I)=8, etc. 

Note: If C(I)^ 0 the program 

sets it to 2 . The dependent 
variable has its C(I) set to 
1 by the program. 


SUBROUTINE TRANSG (X, NOV, NVA) 

DIMENSION X(l) 

Statements defining added variables or modifications of 
original variables 

RETURN 

END 
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V. OUTPUTS OF THE PROGRAM 


Output from the program consists of a title page with 
input data information (Figure 5.1), one or more pages of 
the first twenty time history points (Figure 5.2), one or 
more pages showing the mean and standard deviation of the 
variables (Figure 5.3), one or more pages showing the co- 
variance matrix (Figure 5.4), one or more pages of the cor- 
relation matrix (Figure 5.5), and one or more pages of a 
summary for each iteration done by the program (Figure 5.6). 

In addition to the standard output above, there are three 
optional output capabilities controlled by variables on card 
type 12: 

1) A summary table indicating what variables were 
entered or removed from the regression equation 
during the run. (Figure 5.7). 

2) A list of residuals consisting of actual and com- 
puted values for the dependent variable, the 
difference between the actual and computed values 
of the dependent variable, the relative error of 
the difference, and the first five variables 
included in the regression equation. (Figure 5.8). 

3) A printer plot of actual and calculated time 
histories for the dependent variable (Figure 5.9). 
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VI. FILE DESCRIPTION 


Besides the card reader and printer, the following 
files are used by the program: 

TAPE2 contains time histories of the original 
and also the transgenerated variables. 

It is needed to calculate the residuals 
of the variables being regressed. 

TAPES contains the summary table for the current 
regression run. 

TAPE4 contains the time history input data. 

TAPE 9 contains the mean and covariance matrix 
for all variables, both original and 
transgenerated. 
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Figure 5.2 Time History Table 
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Figure 5.3 Means and Standard Deviations 
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Figure 5.4 Covariance Matrix 
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Figure 5.5 Correlation Matrix 
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Figure 5.7 Variables Entered or Removed 


LIST 

CASf 

UUySfK 


OF RESIDUALS 


V 

M U . 






-tmoiiaa iL:!i. 


JiL.lL. 




JLilii- 


>IU4L 


“lilln-oa 

•S^SiC- 0 ^ 
-fC? U-f i 


• jLtOi 



— Tf i*yji_*y5 ! 

*n ••111 

:Oi r4iyi*:.ryi_ 

-02 -.jaw 





-. 7595F-03 

-.16*;9(-02 

• *.)y3i -33 

-U21L 


-•l;i6i-(2 

-.2U7t-02 

•4613C-03 

-.3103 

13 

-.15501-02 

-,2422* -02 



14 

-.2:4<';-i.2 

-•>5 721 -i'2 

• :2J4i.-7) 

J6*: * 

>L- 

-•2664f-02 

-£26/31-02 

ie6r-i,-)5 

-£ u* 7 ai- 

i6 

-r46d5c-02 

-.3^67-02 

-.6 24 17-33“ 

• 25i> 


74 -Wbsa-.^OZ 

"26 -;63J2l-02 

2 ? T» 


_j 7 3701 - 
• l Ooii *-i i 




NJ 

29 

-•401 jJ-02 

•71111-02 

Ov 

10 

-.7046i;-02 

•9054J-02 


31 

-.65971-62 

•6l33fc-02 


ll 

-.l54Vf-01 

•9647C-03 


33 

-.24151-7 1 

-•Icc-iJ-.M 


34* 

-. 3l22f-6l 

-•lV69l-0i 


35 

-.3455^-01 

-.52*i7l-01 


36 

-•356U-01 

-•6234C-01 


• • « • f. •> 



• 16tU-31 2. 2nd 

ryj liUl, 

►n 1-C62 

“ *2 £ 26.. It 

-52 *'-.ujr 

lil 



46 

47 
46 
49 

" 56 “ 


• 3I04>0^ 

rjs74i*-.M 

. 2137 ^ 0 ! 

• 2l44t-0r 

• 2<i3PE-ai 

Talt’-ii 

• 3j9Vl-f • 

;i5i5f-tL' 

• 3(^*0F^^l_ 

ri 7 t 4 ;* 0 | 

• 37i^f:*0| 
-.3t^0i-01 
-.1263t-0| 
- 72966:^1 


-•^ 4 i:?r-rr 

_^,^72 4^-( I 
-^ 4 nu- 6 i 
-, ^ C|0|-01 
-•40041 ‘ 

— • <*_42 1 

■-i 404 uZ 

-♦ 3***>il 

-ni 64 r 
^ziUddi 


n -•7509 

n 

U -K192 
J1 il-2H 


t?ii7?.ryi .... mUtrOa. 

>J2 •IOU9F-Ot •i727Je-0| •2Jid£-02 •1102 

•.-yi f .ikil .iiaUT))2L #4062. 

-Jl .l7-*»'>i. •i4w9 *72iili-l2 .4Vrf9 

iJiiJ ^Uhilzll .aliiirJZ 


-01 -.|j5i£-0l 

-n -•914C*f-0l 

?nir:.-n u .;;2 

• llt«jL-^i -.154b 

T. 2 QIL 



is 

-•164df-H 

-.IWU-Dl 

-.52l2t-02 

. 3423 

19 

-tl?rs!;rLl_ 

-.72v3^-;*2 


• 466 1 

iQ 

-.i525i-6r^ 

-•If vil -31 

-.47iil-02 

• 3o;4 

21 

-.I372f-0l 

-•62J(** -02 

-,7432f-02 

• iU;* 

22 

-;i2i6£-01 

-•147bi-u2' 

-. tUbi'.-3l 

il754 

23 

-.1052f-01 

•2C62I-02 

-• 125*t-0l 



.2241 

ri? 37 l 

-.2U6 

. t IW? ir*il r. mtit. 

. J/yii:-ll -.1975 

>f.)lQ?Lryi r.UiU 

•i/bdu-Ol -•7<d7L-Jl 

_ j W66t-0| *9t»3^i 

•i2v.c-6l .3260 

^j4?!s7£ry| 

• ll:i2i.-0l .ladQ 

_xlilri^:r \t-Zt1- 

• d619<.-u2 K679 

2,m_ 

-.iioOl-02 2,575 


tio;> 

• 44)3 

• i J60 

• I23C4 

• 2710 

• •3v6i 

• 1 Hi ^ 

• 2234 

• iiw 

• laud 

12Q5_ 


•9720C-02 

_.4QJJ;jrOl_ 

•I103C-31 

^iUoqiz 9 \. 

•1)145-01 


• 71 «U 

-.7462 

• 6 ) cb 

..i9V(»5 

• udui 

W40i 


•l390i-01 «dl04 

_.i720;.TiL .7c70_ 

• lt4.i-01 •7146 

.tX)))s:Oi .654tt. 

•i9l7c-Oi •5676 

^^AUlizn .,26».. 

•ilfti-ll -•2C64 


• t*db3..-01 


_.il554i-il. 

• 1477:-01 

• I1474C-01 

•1)57;-U1 

.-•W70ir0l 

• I W2i-Ul 


. r.4)lO... 

-•66C6 
-r.570Q_ 
-• i4oO 

-1.20b 

-1. 3od 



-. I jCu6-:ii 

3.013 

_,tblb';OJ 

i>)292:*y 1 

-1.15) 

-.2iwt-n 

1.4f4 

• lv71 

. 106 Ji-Ol 

-. d4 95 


- i.or4. 

.1712 


“t 77/1. 

-.b34r:-Ol 

4.275 

• b76U-n 

• bb*«2:*i2 

-.5352 

“1 Jv 7oer01 

_5i^44 


• 643) s -02 

-• 3077 

-.1267 

4.664 

-•51<.it-01 

.33/9c-j2 

-•4C.92 

-.147? 

4.625 


-r454Vv-.}2 

-• dclb 

-.1177 

4.441 

-.446JE-U1 

-.13465-31 

-1.015 

-.1946 

4.151 

-* ?4 w 1 C-5| 


**i t204 

-• 1 nt 

3.746 

-.1425t-Ol 

-.4025^-01 

-..34) 

_::.16V3 

.]f401 .. 




—.1436 

3.002 

.29251-01 

-.2)11^-01 

-l.titf 

^iU?2 

2£5(>9 

-.24V4l:-Oi 

-•2792i-3l 

— • 70S4 

-• 1376 

2.116* 

-. 115u 

-.2/766-01 

-.7d/9 

-.9l/3:-01 

l.55n 

-.1 767 



Jl 

i.’Ia 

-tiiiu 

-• 4 li 4£-4| 

-• U4'|7 

-.99125-01 

t1Vi2 

-•290) 

*iI|74(-^1 

- » *4 4 6 

-•10^3“ 

.4639 

-. 3569 

-.43301-01 

-.55/2 




-.•195o __ 

^t4740«r0l 

_r% 534^ 

-•1342 

-.71632-01 

-.4334 

-.60744-01 

-.5179 


-.2175 

-.4541 

> > J /*: 1 


-.U71 

-.424^> 


t 1 & -I n I • ^ ^ 1 1 

-•5427C-OI 

-.4464 

-•1)35 

-.5777 



^•4497 

-;i96r" 

-.7330 

-. 5355 

-•5l39i-0l 

-.4646 


Figure 5.8 Residuals 
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Figure 5.9 Plot of Regression Run 





